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Abstract 

A self-consistent extended Einstein-Maxwell model for relativistic non-stationary pola- 
rizable-magnetizable anisotropic media is presented. Based on the analogy with relativis- 
tic extended irreversible (transient) thermodynamics, the extended constitutive equations 
for the electrodynamics of continua are formulated phenomenologically, the convective 
derivatives of the first, second, etc. orders being taken into account. The master equa- 
tions for the gravity field contain a modified effective (symmetric) stress-energy tensor of 
the electromagnetic field in a material medium, the use of this tensor being motivated 
both by historical analogies and direct variational procedure. By way of example we con- 
sider the exact solution of the extended Einstein-Maxwell model, describing the isotropic 
cosmological model with hidden non-vanishing electromagnetic field, electric polarization 
and magnetization. 
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1 Introduction 



Non-stationarity of electromagnetically active media is known to lead to complicated non- 
equilibrium phenomena, and one of the cooperative degrees of freedom, excited by the medium 
dynamics, relates to the phenomena of polarization and magnetization [TJ-[6]. The Universe is 
a non-stationary system that contains numerous electrodynamic subsystems of different scales 
[7J [H [9]. It is natural to expect that the cosmological non-steady background stimulates var- 
ious evolutionary processes, involving the phenomena of polarization and magnetization dy- 
namics. This circumstance attracts our attention to the exact non- stationary solutions to the 
self-consistent Einstein-Maxwell equations, describing the evolution of gravitating polarizable- 
magnetizable media. Following basic Einstein's ideas, the stress and energy of polarization and 
magnetization themselves act as sources of the gravity field. These sources can not be negligi- 
ble in comparison with the contribution of pure electromagnetic field: for instance, the plasma 
susceptibility parameters can be at least of the same order than one (or greater), i.e., the contri- 
butions of the electric polarization and pure electric field to the electric induction of the plasma 
are comparable. Thus, in order to formulate a self-consistent Einstein- Maxwell model one needs 
to find the adequate energy-momentum tensor of the electromagnetically active medium taking 
into account the polarization and magnetization of the medium. 

In many cosmological models the total stress-energy tensor is presented as a sum of the 
stress-energy of the fluid (perfect or viscous) and of the energy-momentum tensor of pure 
electromagnetic field [TUl ITT] . To take into account the interacting (cross) terms, which in- 
clude the polarization and magnetization of the medium, one needs to overcome the so-called 
"Minkowski- Abraham controversy" , which followed the famous papers of Minkowski [12] , Ein- 
stein and Laub [13] and Abraham [T3]. In the second half of the last century the problem of 
representation of the stress-energy tensor for the electromagnetically active media revived. The 
modified versions of the stress-energy tensor, related to the specific properties of the so-called 
ponderomotive force, have been presented and motivated by de Groot and Mazur [15], Grot 
[T6] . Israel [17], Maugin [18] and others (see, Refs. [19]- [23], [5], [6] for reviews). Some new 
aspects of the interrelation between the Minkowski and the Abraham versions of the electro- 
magnetic stress-energy tensors and their modifications, can be found in [24|-[31j. However, 
the solution of the 'Abraham-Minkowski controversy" might lay in a modification of termi- 
nology. Particularly, Gordon [29J introduced the label "pseudomomentum" in application to 
the Minkowski momentum. Nelson motivated the use of the term "wave momentum" for the 
sum of the Abraham momentum and of the Minkowski one [30] . Garrison and Chiao [31] used 
the terms "canonical and kinetic form" of electromagnetic momentum. It is clear now, that 
there are two different aspects in the problem under discussion. The first is connected with the 
correct definition of the flux four-vector of the electromagnetic field in the material medium 
as a part of the so-called electromagnetic energy - momentum tensor [18J, which appears in 
balance equations. The structure of this tensor can be verified in the laboratory, and a few 
experiments have been proposed for this purpose (see, e.g., the review [23])- The second aspect 
is connected with the correct construction of the so-called effective stress-energy tensor of the 
electromagnetic field as an adequate part of the total stress-energy tensor, which appears as 
the source of the gravity field in the self-consistent Einstein-Maxwell model. It is interesting 
to clarify the structure of this tensor for the cosmological applications and to check it using 
the observational data. Here we focus on the problem of representation of the second quantity, 
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namely, on the effective stress-energy tensor of the electromagnetic field in a polarizable and 
magnetizable medium. 

The paper is organized as follows. In the Section 2 we discuss the reconstruction of the 
electromagnetic energy-momentum tensor on the base of balance equations, which are the con- 
sequences of the Maxwell equations; as well, we concern some historical aspects, which are 
necessary for further consideration. In the Section 3 we introduce a new version of the effec- 
tive stress-energy tensor. For this purpose, using Lagrange formalism, we derive directly this 
tensor for three well-known cases: vacuum, spatially isotropic medium and uniaxial anisotropic 
medium. Then, based on the obtained exact formulas, we propose an ansatz about the structure 
of the effective stress-energy tensor for general case and compare it with well-known tensors 
of this type. In the Section 4 we discuss the stationary and non-stationary phenomenological 
constitutive equations, linking the polarization - magnetization tensor and its derivatives with 
the Maxwell tensor and its derivatives. Thus, we introduce the so-called extended constitutive 
equations for relativistic electrodynamics of continuous media, and in this sense we get an ex- 
tended Einstein-Maxwell model. In the Section 5 we discuss the example of the exact solution 
of the Einstein-Maxwell equations, based on the proposed effective stress-energy tensor of the 
electromagnetic field. To demonstrate the novelty of our approach, we discuss a particular 
example, when electric and magnetic fields, polarization and magnetization, as well as electric 
and magnetic inductions in the medium are non- vanishing, nevertheless, their cooperative con- 
tribution to the total stress-energy tensor is equal to zero. In this exact model the Einstein 
equations coincide with the ones from Friedmann-Lemaitre- Robertson- Walker (FLRW) model, 
the electromagnetic field and electromagnetic induction being hidden from the point of view of 
space-time evolution. Section 6 summarises our findings. 



2 Electromagnetic energy- momentum tensor 
2.1 Balance equations 

The standard phenomenological way to introduce the electromagnetic energy - momentum 
tensor in a material medium is connected with balance equations (see, e.g., [32] for details). 
The balance equations can be derived using the Maxwell equations 

•> 47T ■ 

V k H lk = P, (1) 

c 

ViF kl + V t F lk + V k F u = , (2) 

where F tk is the Maxwell tensor and H lk is an induction tensor [TJ [2j [33] . The convolution of 
equation flTJ with F\ 

F\S7 k H ik = FF\ (3) 

c 

can be transformed, using (j2J), into equations, which have an explicit divergence form 

V fe T M = F l . (4) 

The structure of (j4]) allows one to indicate the pair T kl and F l as a conjugated one, T kl being 
an electromagnetic energy-momentum tensor, and F l being a ponderomotive force. The choice 
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for T kl and F l is not unique. To fix this pair one should use some supplementary motivation 
(microscopic or phenomeno logical) . Consider now six well-known instances. 

2.1.1 Vacuum case 

In vacuum (P = 0, H lk = F tk ) equation takes the conservation law form 



V^j = 



where 



7'kl — n kl rp T?ran jpkm jpl 
(0) — H9 fmn-F — r f 



(5) 
(6) 



__ 1 
: (o) = i' 

is the stress-energy tensor of the electromagnetic field. It is symmetric, traceless, conserved 
and does not depend explicitly on the velocity four-vector U l (of an observer or of the medium 
as a whole). The ponderomotive force in vacuum vanishes. 



2.1.2 Minkowski version 

By rearranging equation as 
1 



V fc 



— n TJ rprnn jrKm rp 



c 4 



F mn V H 



l zjran 



I T?mn 



TJ V7t p 



(7) 



we obtain the conjugated pair T and F 



rykl 

1 (Minkowski) 



„kl tj T?mn Tjkm t?1 
tt mn b - ti t m , 



(8) 



r (Minkowski) — H r ' ^ 
^ik 



F mri V l M mn - M mn V l F mn 



(9) 



Here M = H — F is the polarization - magnetization tensor of the material medium. The 
tensor ^Minkowski) ^ s traceless, but not symmetric. It does not include explicitly the velocity 
four-vector of the material medium. 



2.1.3 Modified Minkowski tensor 

Grot and Eringen [31], Israel [17], Maugin [18] , proposed a modified version of T kl and F l in the 
balance equation (j3J). These authors discussed the following electromagnetic energy-momentum 
tensor and ponderomotive force 



rrkl 

i (modifl) 



-a kl F F r 

A ran- 1 - 



jjkmpl 



Air 

F l = —InF™ 1 



_ M mn v lp 



This tensor, whose trace does not vanish, T, 



kl 

(modifl) 



c 2 
differs from (181) in the first term. 



(10) 
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2.1.4 Version of Hehl and Obukhov 



Hehl and Obukhov motivated in [28] the following choice for T kl and F l 

^ldif2) = \g kl F mn F mn - F km F l m , F l = ^I m F ml + F l m V k M km . (11) 

The structure of this stress-energy tensor coincides with that of the vacuum, thus, it is sym- 
metric, traceless and does not depend on U l . Nevertheless, in contrast to the vacuum case, 
it is not a conserved quantity, and the ponderomotive force is linear in the divergence of the 
polarization- magnetization tensor. 



2.1.5 Abraham's version 

Abraham [14] proposed to use the symmetric electromagnetic energy-momentum tensor, which 
depends explicitly on the U l , velocity four- vector of the medium as whole. When the medium 
is spatially isotropic and homogeneous, the corresponding tensor ^Abraham) rea ds 

^Abraham) = ^(Minkowski) + ~ ty^U , (12) 

where n is a refractive index of the medium and the vector fl 1 is 

Q* = Ui(H il U m + H^U 1 + H mi U l )F ms U s . (13) 

Here and below we shall use the normalization U k Uk = 1. Since Q l Ui = the tensor ^Abraham) 
is traceless. The corresponding ponderomotive force can be obtained from (j4"l). 



2.1.6 Version of de Groot and Suttorp 

De Groot and Suttorp [36J , based on a microscopic motivation, proposed the following electro- 
magnetic energy-momentum tensor 

^Idifs) = ^Idifi) - U l U m (F kn M nm - M kn F nm ) + U k U l U m U n F ms M sn • (14) 

This tensor also depends explicitly on the velocity four-vector. The corresponding ponderomo- 
tive force can be obtained from (j3J). 

Remark on the microscopic and macroscopic electrodynamics. 

As it was emphasized, e.g., in [361 El], the problem of the choice of the electromagnetic 
energy-momentum tensor T kl is connected with the basic microscopic model, as well as, with 
the averaging procedure of the microscopic Maxwell equations (the discussion about averaging 
procedure see, e.g., in |35j). Indeed, let the microscopic electromagnetic field f ik can be repre- 
sented as a sum of a mean field F^ and a fluctuation terms with vanishing average value 
(£ifc) = 0. Then, one obtains, that 

i 

Tik / ik f rmn rim ek \ rpik , l~.ik\ /ir\ 
= JmnJ -J I ml — J (0) + \ T /) 

where TAft is given by ([6]). The second term 

/_ife\ / „ifet tmn s-im/-k \ / t cz\ 

\T ) = {^9 Un4 -4 Cm) ( 16 ) 

is very sensible to the averaging procedure and essentially depends on the microscopic model 
of electromagnetic interactions in the medium. 
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2.2 "DEHB" - representation of the energy- momentum tensor 

The tensor T lk can be rewritten in terms of four- vectors D l , E\ H l and B\ which play signif- 
icant role in the covariant electrodynamics of continuous media [H]. The definitions of these 
four-vectors are well-known [37] 

D l = H lk U k , Hi = H* k U k , E^F^Uk, B l = F* k U k . (17) 

D\ E l , H l and B l are orthogonal to the U\ a four- vector of macroscopic velocity of the medium. 
In its turn, F lk and H lk can be represented as 

pik = E ijjk _ E k v i _ v ikj B ^ ^ H ik = D i uk _ D k v i _ v ikj H , ^ ^ 

where 

rpikjs 

v ikj = e ikjs Us ^ e ikjs = __ _ ( 1Q ) 

e ikjs - g ^ e Levi-Civita tensor and the term E lk ^ s is the completely skew-symmetric Levi-Civita 
symbol with E 0123 = 1. This tensor provides the dualization procedure: F* k = \^ikmnF mn ■ By 
means of (JTSJ the tensors ^ inkowski) , 7^ braham) , 7^ 0<m} , T™ odi{2) , 7g? odif3 ) can be represented 
in terms of four- vectors D l , E\ H l and B % . For instance, the Minkowski tensor has the form 

Minkowski) = - U p U^ (D m E m + H m B m ) - 

(D p E q +B p H q ) -U v rf mn E m E n -U q rf mn D m B n . (20) 
Taking into account the standard decomposition of this tensor 

^(Minkowski) = W( cm )U l U k + U l I k ^ + U Ify + ^(Minkowski) J (21) 

one can conclude that the energy density scalar Wr em ), the first and second flux four- vectors 
P {1) , P {2) and the stress tensor P ( M inkowski ) read ' respectively, 

W (em) = U p T$ mkowski) U q = \ (D m E m +H m B m ) , (22) 

^(i) = ^p^Minkowski)^? = ~V k mn D m B n , (23) 
I{ 2) ee A;^ inkowski) [/ 9 = -v l mn E m H n , (24) 

^Minkowski) = A;^ inkowski) Aj = ~ A* k (D m E m + H m B m ) - + B*H k ) . (25) 

Note that the Minkowski ([8]) and Abraham ffl2|) versions of the energy-momentum tensor share 
the same quantity Wf em \. For the Abraham (symmetric) version of the electromagnetic energy- 
momentum tensor, Ihs and /Ln coincide and 

= 4) = "^n™" • (26) 



6 



The stress tensor ^(Abraham) coincides with the symmetrized one, which was obtained by Min- 
kowski [23J. In terms of three- vectors E, B, D and H the flux three vectors proposed by 
Minkowski, read, respectively, 

I (1) = [D,B], I {2) = [E,H), (27) 

where [D, B] denotes the vectorial product of the three- vectors D and B. In spatially isotropic 
medium one has D = eE and B = fiH, where e and /i are the scalars of electric and magnetic 
permeability, respectively Thus, using the standard definition for the Poynting flux three- 
vector, S(poynting) = [E, H], one can write for the corresponding three- vectors of the momentum 
of the electromagnetic field 

^(Abraham) = ^(Poynting) > ^(Minkowski) = ^(Poynting) • (28) 



3 Effective stress-energy tensor 

3.1 Lagrange formalism 

The Einstein field equations 

Rik _ }_ g ik R = Ag ik + KT k^ j (29) 

must have on their right-hand side the so-called total stress-energy tensor TA„ ta n, which must 
be symmetric by definition and divergence-free due to the Bianchi identities |38j, i.e., 

rpik _ rphi v-7 rpik f\ ( OC\\ 

1 (total) — 1 (total) > * fc (total) ~ U • i 6U ) 

R tk is the Ricci tensor, R is the Ricci scalar, associated with metric and A is the cosmological 
constant. Following the standard variation procedure one can define Tu ota] ^ as 

T ^'> s -T^i ^ ' (31) 

where the scalar £ denotes the Lagrangian of the whole system, and includes the terms related 
to the electromagnetic field, the polarization and magnetization of the medium. The main 
problem is how to separate the contribution of the pure electromagnetic field, the contribu- 
tion of the polarization and the magnetization and the contribution of the pure matter. This 
problem seems analogous to the problem of separation of pure gravitational energy-momentum 
and the energy-momentum of the medium, which is characterized by the gravitational self- 
interaction. One can extract from the total stress-energy tensor T^ otal ^ the electromagnetic 
energy-momentum tensor T lk (see, previous Section), which, in general, is not necessarily sym- 
metric and traceless. On other hand, based on the variation procedure with respect to metric, 
we can draw from T^ otg ^ the so-called effective stress- energy tensor of the electromagnetic field 
which is symmetric and traceless by definition. Below we will distinguish between T %k 
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and Tu S y In order to motivate our ansatz about the effective stress-energy tensor, let us, first, 
consider the variation procedure of its derivation for the simplest action functional 

S[F mn ,g pq \ = J d 4 xy/^g j — h L( ma tter) + -jC lkmn F ik F mr }^ . (32) 

Here C tkmn is the linear response tensor, which describes the influence of matter to the electro- 
magnetic field. This tensor has the following symmetries 

^jikmn ^jkimn ^jiknm ^jmnik (33) 

Variation of S[F mn ,g pq ] with respect to the four-vector electromagnetic potential Ai gives the 
Maxwell equations 

V k H ik = 0, H ik = C ikmn F mn , (34) 

where H lk is the induction tensor and the current four- vector P is absent. The variation with 
respect to metric tensor g pq yields the Einstein equations fl29l) with explicit decomposition 

1 (total) ~~ 1 (matter) 1 (eff) ' J 

As usual, the symmetric stress - energy tensor of the material medium Tf^ atter -) reads 

T {Ltter)=WU i U k + q l U k + q k U l - PA tk + Tl ik , (36) 

where W is an energy density scalar of the matter, IP is a macroscopic velocity four-vector of 
the medium as whole, q % is a heat-flux four- vector, P is the Pascal pressure, A tk = g tk — U l U k 
is a projector and IP fc is an anisotropic pressure tensor. The form of the electromagnetic part 
of the total stress-energy tensor, TjPjjL, depends on the suggestions about a structure of C lkmn 
tensor. When C tkmn incorporates the metric g pq only, the effective stress-energy tensor takes 
the form 

rppq _ 1 jpqr^ikmn p p _ jspqikmn p p /q7\ 
(eff) — r ik r mn g r ik r mn i \° ' ) 



where the tensor K pq mn is a formal variation derivative 

j^pqikmn _ ^jikmn (38) 

5g P q 

is, by definition, a symmetric tensor, whose trace vanishes when 

In this paper we assume, that the tensor C lkmn contains the metric only. When C lkmn contains 
the Riemann tensor, the Ricci tensor and the Ricci scalar, the corresponding effective stress- 
energy tensor includes the covariant derivatives of the Maxwell tensor up to the second order, 
and we deal with the so-called non-minimal Einstein - Maxwell theory [39l HQ]. When C lkmn 
includes the covariant derivative of the velocity four- vector Vjt/fc, the effective stress-energy 
tensor involves first covariant derivative of F mn and we deal with dynamo-optical effects [41] . 
But it is worth stressing once again that here we restrict our-selves by the first case only, 
namely, when C %kmn = C tkmn [g pq \. 
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3.2 Structure of C ikjl 

When C 3 = C 3 , a standard decomposition of the C 3 tensor in terms of dielectric per- 
meability tensor e tk , magnetic impermeability (ji~ 1 )ik tensor and magneto-electric tensor v k is 
permissible 

^jikvan 
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^.imjjkjjn ^injjkjjm _|_ ^.knjjijjm ^kmjjijjn 



\rf kl {^r 1 )i s rf nns -]- [r 1 ik \U m v l n -U n v l m )+rf mn (U l v l k -U k v l l )\ , (40) 
2 2 ~ 

where 

£im = 2C ikmn UkUn ? = Arj pik C ikmn r, mnq , v v m = V pikC ikmn U n . (41) 

The tensors and {[i~ l )ik are symmetric, but v k is, in general, non-symmetric. These three 
tensors are orthogonal to U l , 

6^ = 0, {^ l ) ik U k = 0, v l k TJ i = Q = v l h U h . (42) 

The case C 3 ^ (jjhk j g d escr ib ec i j n detail by Hehl and Obukhov [33] • The decomposition 
fT40|) leads to the well-known formulas for the four- vectors of electric induction D l , magnetic 
field Hi, electric field E l and magnetic induction Bf. 

D i = e lk E k - v k i B k , E i = v l k E k + { l i~ 1 ) lk B k . (43) 

The material tensors e lk , (/x -1 )^, v x k and C lkmn can be decomposed using the standard tetrad 
representation 

ci\ii—%n — JSf»l . . J^in c(ai)(a2)--(a n ) ('44") 

(01) (02) (dn) ' v ' 

Here the symbol .X? \ denotes the set of the four tetrad vectors, whose index (a) runs over 
(0), (1), (2), (3), and XL\ = U l . These four four- vectors are assumed to satisfy the orthogonality 
- normalization rules 

9ikX\ a ) X {b) = V(a)(b) , (45) 

V (a)ib) XU X (t» = 9 M > (46) 
where T]< a ){b) denotes the Minkowski matrix, diagonal (1,-1,-1,-1). Since the tetrad four- 
vectors are linked by the relation containing the metric, for further consideration we have to 

define the formula for the variation , . 

3.3 Variation of the tetrad vectors 

Varying the relations (j4"5]) with respect to the metric, we obtain 

X k (b)5X k a) + X k ( a) 5X k b) = -X^X^dgn, . (47) 
The variation of (T46l) yields 

Sg» = J*)® [x« b) 6Xf a) + Xf a) 5Xfa] . (48) 



The relations (|47p and (j4"81) are equivalent, since 



9 ik 9kj = &j => 5g ik = -5g pq g pi g qk . (49) 

The variation of arbitrary origin SX 1 ^ (not necessarily caused by the metric variation) can be 
represented as a linear combination of the tetrad four-vectors: 

8X\ a) = X{ c) Y { / ] . (50) 

The tetrad tensor K% is not generally symmetric. Using the convolution of (148|) with tetrad 
vectors, we obtain 

Y (a)(b) + Y (b)(a) = 5 gPQX^X^ , (51) 

where we use the standard rules for the indices, e.g., X^' = rj^ h ' g qm X^y Consequently, the 
symmetric part of the quantity Y^ b \ indicated as Z^ h \ can be readily found: 

Z (a)(b) = ^gPQX^X^ , (52) 

and the law fl50|) reads now 

5X\ a) = hf* [X p{a) 5 q + X q(a) 6 p ] + X\ C) Z$ . (53) 



Here Z'$ is the skew-symmetric part of Y^j, i.e., 2i£( )( c )=Y( )( c )— Y( c )( a ). Therefore, the vari- 
ation of the metric produces the variation of the tetrad, described by (153]) with vanishing 
skew-symmetric part Zu\i c \. Thus, one finally has 

ML 1 r i 

" -- [X pia) 6l + X qia) 6;] , (54) 



Sg 

and we can use this formula for the variation of the four- velocity vector U l and for the variation 
of the space-like vectors XJ a ^ (a = 1, 2, 3). 

Note that in [12] the author, considering the stress-energy tensor for the spinor field, has 
used the formula for the variation of tetrad, which can be easily transformed into (I53p . The 
authors of [23] - [27] use a different formula for the variation of the four-velocity vector: 

/ fir^ 1 

5U\s) = 5 - = --U^g pq UW . (55) 

This formula is derived on the base of kinematic representation of the four- velocity U l (s) as a 
tangent vector for the observer world-line. Such a representation is adequate to the Lagrangian 
variation procedure with respect to coordinates. However, it is not possible to use the same 
method for the calculation of the variation of the space-like vectors XJ a y 

Our approach is based on the consideration of the four vector fields XI J x), subjected 
to the orthogonality and normalization conditions. This representation is appropriate for the 
procedure of variation of the Lagrangian fl32l with respect to the Ai(x) and g pq (x) fields, as 
well as with respect to spinor field, as it follows from 
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3.4 Three standard examples 

3.4.1 Pure vacuum 

The tensor C lkmn of vacuum must be constructed from the metric only, 

C ikmn = 1 (gimghn _ g in g km ^ . (56) 

Then 

jy ikmn ^ ( cm kn i zk <rn im ci cn„km srfc <rm „ik\ /"KT*\ 

=2VV^ + d P d g9 -o p d q g -d p b g g ) , (57) 

and the direct variation in (138]) gives the standard formula for the symmetric traceless stress- 
energy tensor of the electromagnetic field in vacuum ([6]). 

3.4.2 Spatially isotropic medium 

By contrast to the vacuum case the Lagrangian of electromagnetic field in the spatially isotropic 
medium should contain one supplementary vector, namely, the velocity four-vector of the 
medium, U\ The linear response tensor has the following form [T81 15]: 



^jikmn ^ 

~2/i 



(g im g kn -g in g km ) +(e/i-l) (g im U k U n -g in U k U m +g kn U i U m -g km U l U n )] , (58) 



where e and fi are the dielectric and magnetic permeability scalars, respectively. The e lk , 
(A* -1 )?? an d v\ tensors, entering the linear response tensor (see, (|4"T1) ) read 



£ 



1 



"'' - sA ik , (^) M = -A M > "\ = 0. (59) 



Since the velocity four- vector is normalized as UkU k =l, we should use the variation of £/* with 
respect to g pq in order to obtain the K pqikmn tensor (1381) . To do this we employ the formula 

fl7< = Ug™ (u p 5 q + U q 5 p ) (60) 

as a particular case of (I54|) with (a) = (0). A straightforward calculation shows that the 
corresponding effective stress-energy tensor can be written as 

Tgjj = ^C^F™^ - i(C fcm ^F' m + C" m ^F fc m )F pg . (61) 
3.4.3 Medium with uniaxial symmetry 

The tensor C tkmn for the uniaxial symmetry contains not only the velocity four- vector U\ 
but one supplementary space-like four-vector, (say, X 1 ) as well. This vector is normalized 
according to gikX l X k = —1, and orthogonal to U l , i.e., gikX l U k = 0. Thus, to calculate the 
tensor x pgikmn we must find the variation of X 1 with respect to the metric g pq in addition to 
the variation 5U l . For this purpose we use the formula 

6X* = -5g™ (X p 8\ + X q 6 p ) (62) 
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as a particular case of fl54|) . To modify the C lkmn tensor for the uniaxial case we follow the 
procedure described in [13]. First, we modify (1591) by 

s ik = e (A ife + £X*X k ) , ^-% q = - (A pq + (X p X q ) , (63) 

and consider the medium without magnetoelectric cross-terms, i.e., require the following rela- 
tion to hold v tk = 0. Here £ is a coefficient of anisotropy of the dielectric permeability. In the 
uniaxial case it is unique, and can be defined, for instance, as £ = 3 — s\je. This parameter 
vanishes when medium is spatially isotropic. ( is a corresponding coefficient of anisotropy of 
the magnetic permeability. Then, we use (I63p in (l4"Uj) . and after a long but otherwise straight- 
forward calculation we obtain the expressions for j{PQ tkmn anc l Tf£ S y The result was to be 
expected: we recover the formula fIBTl) for the TfJL. 



3.4.4 General case 

In general, the tensor of linear response can be represented as 

C ikmn = xfaXfaXfoXfaCWMW . (64) 

Our ansatz is that the tetrad quantity C^ b ^°^ does not depend on the metric, and the 
variation with respect to g pq reduces to the variation of the tetrad four- vectors only, i.e., to the 
formula (1541) . Under such an assumption using f[5"4"|) and (Hoi) we obtain again the expression 

(EU). 



3.5 Ansatz 

Based on the direct derivation of the tensor T k ^ together with the Lagrangian (1321) . and taking 
into account the common structure (1611) for all three well-known models with the constitutive 
equations H lk = C lkmn F mn , in which C lkmn contains the metric only, we propose to consider 
the following effective stress-energy tensor of the electromagnetic field with general constitutive 
equations 



1 1 



Tfa = -g kl H mn F mn --(H km F t m +H lm F k J . (65) 

If we use the effective stress-energy tensor T k ^ as an electromagnetic energy - momentum 
tensor T kl in the relation (J3J), the corresponding ponderomotive force takes the form 

Fr - ^ + \ - M„„V^, + iv. [M*»F lm - ^M lm] . W 

The effective stress-energy tensor ( 165|) coincides with the symmetrized Minkowski electromag- 
netic energy-momentum tensor (JSj) - Thus, it is manifestly symmetric, traceless and does not 
depend explicitly on the choice of U l . 
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V\ k ca ) = AlTtLAj = -A lk (D m E m +H m B m ) — {D i E k +D k E i +H i B k +H k B i ) . (70) 



3.5.1 DEHB - representation of the effective stress-energy tensor 

In terms of four- vectors D l , E\ H l and B % the tensor Tj^L can be represented as follows 

T (es) = i^9 m ~ U p U^ (D m E m + H m B m ) - 

- I (D p E q + D q E p + H p B q + H q B p ) - i (U p r) qrnn + [/y mri ) p m S n + E m H n ) . (67) 

The energy density scalar the flux four vector JL ff -j and the stress tensor T 3 ?^ read, 
respectively, 

W {eS) = U p T pq s) U q = - 1 - (D m E m + H m B m ) , (68) 

P m ee U p T P : s) Al = ApT^U q = -\rf. mn (D m B n + E m H n ) , (69) 

(vid = -y.,,.,,-,-^ - -A lk (D m E m +H m B m ) -- 

Note that the energy density scalar W{es) coincides with W( em ) given by ( |22|) . obtained for the 
Minkowski and the Abraham tensors. The flux four- vector i| eff ) is one half of the sum Ihs and 
IU\ ( |23l) . ( |24l) . The stress tensor Vuq^\ is the symmetrized one, obtained by Minkowski, and 
coincides with the Abraham stress-tensor. In spatially isotropic medium in the three-vector 
symbols we have for our definition of the effective stress-energy tensor 

S(eS) — 2^ £ ^ ^)^(Poynting) = ^(Abraham) = ' ^, ^ (Minkowski) • (71) 

4 Constitutive equations 

The ansatz (joo| concerning the structure of the effective stress-energy tensor of the electro- 
magnetic field in a material medium allows us to consider self-consistently an Einstein-Maxwell 
model for the evolution of non-stationary electromagnetically active media. This model in- 
cludes, first, the Einstein field equations (1291) with (I35I) . (|36|) and (1651) . secondly, the Maxwell 
equations ([I]) and ([2]), thirdly, the so-called constitutive equations, linking the polarization 
- magnetization tensor M lk = H lk — F lk and Maxwell tensor F mn . The specific feature of 
the interaction between a non-stationary material medium and the electromagnetic field is the 
dynamics of polarization and magnetization, and the constitutive equations must reflect the 
existence of an inertia in the electromagnetic response of the medium. In order to motivate our 
ansatz about the constitutive equations, we briefly consider the well-known classical analogs. 

4.1 Classical and relativistic extended thermodynamics as a hint for 
the construction of the covariant extended continuum electro- 
dynamics 

4.1.1 Rheological models 

The well-known Hook law 

a/3 _ na/3jp _ a/3 / *70\ 

a — U e 7P — a (stationary) ' l' Z J 
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describing the linear relation between the stress tensor a af3 and the deformation tensor e 7P , is 
considered to be the simplest stationary linear constitutive law in rheology |4]. Greek indices 
run over 1,2,3 and describe the spatial coordinates. The material tensor C al3lp is symmetric, 
i.e., 

Qftcfyp Qaf3py Q^pafi (73) 

and includes the elastic coefficients of the medium [17]. When the medium is non-stationary, 
the difference a a/S — Stationary) becomes a function of time due to inertia effects. Maxwell's 
viscosity model [H] assumes, that this difference is proportional to the time derivative of the 
stress tensor. In the Kelvin- Voigt model (we follow the terminology of Ref. [45]) the difference 
a ap _ cr?^ atioll s is proportional to the time derivative of the deformation tensor. The Poynting- 
Thomson model (see, [43]) combines Maxwell and Kelvin- Voigt models and is characterized by 
the constitutive law 

U L- C 7p — L yp U T A Q-f- fP 1 " V / 

Here the first term in the right-hand-side corresponds to the Maxwell model and the second 
term relates to the Kelvin- Voigt contribution. The quantities and A a/37P represent the 
tensors of relaxation parameters for stresses and strains, respectively. These three basic models 
provide a rule for the next generalizations of the rheological models (see, e.g., [3]), namely, to 
introduce the successive derivatives of the stress tensor and/or strain tensor (see, e.g., Jeffreys 
model [13]), by the same method, used for the first derivative. 

4.1.2 Heat conduction model 

The Fourier law 

Q= — AVT = ^(stationary) (75) 

is a stationary constitutive law for the heat conduction, connecting the heat flux vector q with 
the spatial gradient of the temperature. Taking into account the inertial properties of the heat 
propagation, Cattaneo [46] considered the difference q — ^(stationary) to be linear in the time 
derivative of the heat flux vector 

q + XVT=-r (q) ^. (76) 

The latter expression, supplemented with the balance equation for the internal energy leads to 
the hyperbolic equation, governing the temperature evolution 

^ t+ ! t =^ 2t - (77) 

which generalizes the standard parabolic equation. By extending of the constitutive law for the 
anisotropic media [47] and using Cattaneo's proposal of the inertia of heat and its extensions 
one obtains the generalized Fourier-Cattaneo law 



f) 8T 
-W*Vf>T - Ql-q? + X a ^p— + ... . (78) 
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4.1.3 Relaxation of the electric polarization and magnetization 



Taking into account the delay in the response of the medium to the applied electromagnetic field, 
one can use the simplest extended constitutive equations containing the first time derivative 

P = X E-T {P) -P, M = X H-T {m) -M. (79) 
These equations are easily transformed into the general relaxation equation 

d 



< = - 

or 



- (i~ ^(stationary)) (80) 



for the polarization-magnetization, see, e.g., [19]. By analogy with rheology, Kluitenberg 
[50~1 |51~1 [52] generalized the equations for polarization and magnetization vectors of classical 
electrodynamics as 

(m)=(i) \ m / H=(i) V a V 

(fe) / a \ ( m ) (s) / a \ M 

M* = c » H f>_ ^ c« M*+ £ ^ (m) - H*. (82) 

(m)=(l) \ t/t '/ M=(l) \ Ul / 

If the medium is anisotropic and there are cross-effects, such as pyro - electricity and pyro - 
magnetism, piezo - electricity and piezo - magnetism [ITJ [53] , magneto - electricity |54| , electro 
- striction and magneto - striction [El [53], etc., then the "source" XpE 13 i n (Ell) and (pH 13 in 
( 182]) must be supplemented by additional terms 



X a pE P + ir a (T - T ) + d%a^ + vfH? + Q% p E^a^ + ... (83) 

QH P => + m a (T-T ) + + v%E^ + R^ p H /3 a' yp + ... (84) 

Here n a and m a are the pyro - electric and pyro-magnetic coefficients, respectively, describing 
the variation of polarization and magnetization produced by deviation of the temperature from 
its equilibrium value, T . The coefficients d% and h% describe the linear piezo - electric and 
piezo - magnetic properties of the medium, respectively. These effects may be induced by 
the stress tensor a 13 " 1 . The tensor z/ a/3 corresponds to magnetoelectric coefficients; if they are 
non-vanishing, the medium transforms electric fields into magnetic fields and vice-versa. The 
electro- and the magneto - striction coefficients Q^ 1P and Rp lp , describe cross-effects involving 
the stress tensor and the electric or magnetic field, respectively. 

When the electro dynamic system as a whole is under motion, the formulas (IBTl) and f[82l) 
have to be supplemented by the terms containing the acceleration vector and the spatial deriva- 
tives of the velocity vector [2]. 



4.1.4 Relativistic fluid in the absence of electromagnetic interactions 

The parabolic equations for the temperature evolution run into conflict with the causality 
principle, since the rate of the temperature propagation is predicted to be unbounded [551 
[55] . As an answer to this challenge, the extended relativistic (or causal, or second order, or 
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transient) thermodynamics was developed. The history, the fundamentals and applications of 
the extended thermodynamics are presented in detail in Refs. [HI [57] and [58]-[ 
The extended covariant definition of heat flux four-vector: 



</' - \TA> ( iv fc T - DU k ) = r {1) AlDq k + 



e + d io 



r (i) \ 

XT 2 J 



35) 



is a relativistic generalization of the definition ( 1761) . Here D = C/ l Vj is the convective derivative, 
Vj is a covariant derivative operator; O = V k U k is the fluid expansion. Both parts of the 
relativistic tensor of non-Pascal pressure 



n 



ik 



ng) + ~A''TJ. U =//,,,! 



obey expressions similar to ( 185 



r (o) 



n-3ce = r (0)J Dn + -^n 



e + ^(iog^ 



36) 



n? 



l (0) - = T (2)^lnK DIi iO) + 

Here the well-known quantity E lfc is introduced by 



(0) 



e + d hog 



T(2) 



S ifc = iA^A^(V m f/ n + V n U m ) - -A lk Q 
2 3 



(89) 



The principal novelty of the relativistic formulas (!85j) - (IHHl) is that they contain the convective 
derivative of the velocity four-vector DU % , the expansion scalar B as a supplementary terms of 
the inertial origin. 



4.2 Remark on covariant electrodynamics of continuous media 

An obvious analogy exists between the constitutive equations in electrodynamics and rheology. 
In this sense, M tk plays in electrodynamics the role of a stress tensor a af3 in rheology, and F lk is 
an analog of the deformation tensor e a p. The main difference is that while electrodynamics deals 
with skew-symmetric quantities, elastodynamics deals with symmetric ones. Following this 
analogy, one can see a correspondence between the four- vector potential A; in electromagnetism 
and the displacement vector V a in classical elastodynamics. Indeed, the second subsystem of 
Maxwell equations in ([2]) leads to the relation 

F ik = diA k - d k Ai , (90) 

which converts it into an identity. Analogously, from the Saint- Venant conditions (equations) 
[53J in classical non-relativistic elastodynamics one obtains that the deformation tensor has to 
be the symmetrized derivative of some three- vector, the displacement vector V a : 

e a = ~(d a V p + dpV a ). (91) 



16 



Likewise, there is an analogy between the constitutive laws in the electrodynamics and elas- 
todynamics. The well-known linear static constitutive equations in electrodynamics says that 
H lk is proportional to F mn 

ti — U fji = H ( stationary ) , (iiZ) 

where the tensor C J is the linear response tensor. In contrast to (1721) with symmetric tensor 
^jq/?7p f|73p relations (I92p contain the tensor C i (1331) with skew-symmetric indices in 
the first and the second pairs. Alternatively, the polarization - magnetization tensor M lk is 
considered proportional to the Maxwell tensor 

M* = X lkmn Fmn = Mftationary) • (93) 

Here ■y/ Lkmn is called the susceptibility tensor, it has the same symmetry of indices, as C %kmn . 
Relation (193]) is, in fact, the analog of the Hook expression (1721) . The covariant phenomenolog- 
ical generalization of the constitutive equations may also be done in terms of polarization and 
magnetization four- vectors P l and M % : 

pi = M ik U k , Mi = M* k U k . (94) 

The equations (193|) yield 

P l = a ik E k - lk \B k , Mi = 1 k E k + [3 lk B k , (95) 

where a tk , (3i k and 7^ can be obtained from (j4"Ij) by the substitution C lkmn ^mn_ R e lativis- 
tic covariant elastodynamics is much more sophisticated (see, e.g., 0, El LTS] and references 
therein), and the described analogy is not so evident. Nevertheless, let us use the main idea of 
generalization of the elastodynamic constitutive equations in order to obtain extended consti- 
tutive equations for the covariant electrodynamics. 

4.3 Phenomenologically extended constitutive equations for rela- 
tivistic electrodynamic systems 

Based on the analogies described in the previous Section we now introduce generalized phe- 
nomenological constitutive equations for non- stationary electromagnetic media. We consider 
three versions of generalization, which, of course, are equivalent, but can be useful for different 
applications. 

4.3.1 The first version 

The first version of the extended electrodynamics of continuous media assumes that the differ- 
ence between M lk and its stationary value x Fjl can be written as 

M ik _ x ^i? 7= r ifc mn Z)M mn +Ai fc mn ^ derivative terms) . (96) 

Here the velocity four-vector enters explicitly the convective derivative D = U k V k only. It is 
clear that the structure of the formula (1961) is analogous to the Poynting-Thomson law (1741) . 
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4.3.2 The second version 

More generally, one can replace the convective derivative, D, by the covariant derivative; then 
we obtain 

^ife _ x ife J 7 i7 ^ 7 =r ifc^^ ^^r, + Aifc ^ v ^ jP ^„ +nlfc? ^ i v ^ ^7^+ ( higher derivative terms). (97) 

In the last case it is convenient to use the standard decomposition 

Vj-17, = UjDUi + Eji + -A^G + u jt . (98) 

The shear tensor E jm was defined in fl89l . the scalar expansion 6 is again VkU k , and the 
vorticity tensor Uji is defined by 

up = ^A™A?(V m U n - V n U m ) . (99) 

Obviously, equation ( 1971) reduces to ( 1961) when 

T% n = U>T% nn , \' h :i„ /'At,. Q^ = WQ tkm . (100) 

In general there exist standard irreducible decompositions of such tensors similar to the de- 
composition of the C lkmn tensor fj4*0|) . The constitutive law (l9Tj) is applicable not only to the 
non- stationary models, but to the non-homogeneous media, as well. In the latter case the 
difference M %k — M^ s fc tationary ^ should be decomposed using the spatial derivatives, A^Vjt, in ad- 
dition to convective derivative D = t/ fc Vfc. Moreover, the presence of higher derivative terms 
in (]9"Tj) gives the possibility to consider cross-terms of the type (ITS]) , described in [ISJ for heat 
conduction. 

4.3.3 The third version 

Consider now the tetrad {!?,Iy}, where (a) = (1), (2), (3). Here X%* = XJ % is the velocity 
four- vector of the medium and the tetrad vectors Xi a \ are connected with three main directions 
in the anisotropic medium. Define the tetrad components of the vectors of polarization and 
magnetization, of the vectors of electric field and magnetic induction, as well the acceleration 
vector as follows: 

P( a ) = PiXl a ) , M( a ) = MjX( Q) , E( Q ) = EiXfa , B{a) = -BjX( a) , (DU)( a ) = X( a) DUi . 

(101) 

Since the vectors P 1 , M\ E l , B l and DU l are orthogonal to the velocity four-vector and the 
vectors X 1 ^ are space-like, we have the inverse decompositions in the form: 

**=-£*W*(«). M*=-Y,M {a) X\ a)) .... (102) 

(a) (a) 

As a main ansatz for the third version we suggest that in the appropriate tetrad the relaxation 
equations take the form 

P («) = a W E W _ 7 _ rff Kp) DP^ + \%DE<fi + $DB<» + r,$(DU)V> + ... , (103) 
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M H = 7 $EW+p$BW-T$ im) DMW + k$DEW + ^DBW+C$(DU)M + ... . (104) 

In such a context we assume, that the Einstein rule for the indices are valid. The term 
T (p)(p) describes the diagonal three-dimensional matrix of the relaxation coefficients for the 
corresponding tetrad components of the polarization vector. Analogously, three independent 
coefficients of the diagonal matrix rfeLj describe three relaxation parameters, which differ in 
general for three tetrad components of the magnetization vector. We assume that different 
components of the polarization and magnetization vectors evolve with its own relaxation time 
parameters. The constitutive equations (11031) and (11041) generalize the formulas (T5TT) and (152]) . 
advocated by Kluitenberg [HIl EE] • 

The first and the third versions of the generalization are equivalent, when the law of the 
tetrad vectors evolution is fixed. To satisfy the conditions of orthogonality-normalization for 
the tetrad vectors, one can use, e.g., the simplest expression 

DX[ a) = n\ Xf Q) , (105) 

where the tensor Qik is skew-symmetric. Then using 

M ik = pijjk _ pkjji _ v ik jMj ^ ( 106 ) 

with ( 11021) and ( 11051) . and expressing the coefficients T tk mn , etc., via the tetrad coefficients 

t 03Xp)' r (/3X m ) "' e ^ c '' one chains the third version (I103I) . (I104I) of the generalized constitutive 
law from the first one (1961) . 



4.3.4 On the electrodynamics of thermo-visco-elastic medium 

When thermo-visco-elastic processes take place, the right-hand-sides of the constitutive equa- 
tions (196]) . (197]) . (I103p . (I104p have to be supplemented by the heat-flux vector q % and its deriva- 
tives D^q\ the non-equilibrium pressure IP fc and its derivatives D^IT k as discussed above. 
In its turn, the constitutive equations for q l and IP fc have to be supplemented by the corre- 
sponding electromagnetic terms M* fc , F mn ,... D^M' lk , D^F mn , etc. This problem, of course, 
requires a special consideration which takes into account the concept of hyperbolicity of master 
equations (see, e.g., [77]). Thus, in general, one obtains a system of coupled extended con- 
stitutive equations for M lh , q l and IP fc , describing the covariant extended electrodynamics of 
continua. 

5 Example of exact solution of non - stationary Einstein 
- Maxwell model 

5.1 Einstein's equations 

We consider the FLRW cosmological model with line element [TUl [38] 

ds 2 = dt 2 - a 2 (t) \(dx 1 ) 2 + (dx 2 ) 2 + (dx 3 ) 2 } . (107) 
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The magnetic field vector, the magnetization vector, the electric field vector and the polar- 
ization vector are pointed along the x 3 axis. For such a "self-parallel" configuration of the 
electromagnetic field the total stress-energy tensor T^ otgl ^, has four non- vanishing components: 

T 0(totid) = W + X > T l(total) = -- P (l)-^ > T 2 2 (total) = _ - P (2)-^ ; T 3(total) = _ - P (3)+^ • ( 108 ) 

Here 

X = \[H l2 F l2 - H 30 F 30 ] (109) 

is the source term related to the electric and magnetic fields, polarization and magnetization. 
We do not specify here the relations between the energy density scalar W and the diagonal 
components of the pressure tensor P^, Pt 2 ) and P( 3 ). They can describe perfect fluid, viscous 
fluid, etc. In a co-moving frame with U % = 5 l , the gravity field equations reduce to the following 
system [38J 

2 2 

= A + k(W + X), 2~+(~) =A- K {Pf 1) +X), (110) 

... 2 2 

2-+(^) =A-k{P w +X), 2-+(-) = A-k(P (3) -X). (Ill) 

The dot denotes the derivative with respect to time. The Einstein equations are self-consistent 
when 

P (1) + X = P (2) + X = P (3) - X = V {isotl) . (112) 

The formulas (11121) guarantee that the global spatial isotropy of the Universe holds. Differen- 
tiating (111 Oft - fillip leads to the conservation law 



W+X+ 



a. 



) [AX+3W+P {1) +P {2) +P {3) ] =0. (113) 

In principle, the matter pressure may be anisotropic, but with P(i) = P( 2 ) = P(3) — 2X, i.e., the 
longitudinal pressure P( 3 ) compensates the influence of the electromagnetic pressure. However, 
the question arises: whether a non-trivial solution exists for which electromagnetic field and 
electromagnetic induction are non-vanishing, but the matter has an isotropic pressure Pm = 
P( 2 ) = P( 3 ) = P. Obviously, such a requirement assumes, that X = 0. For the well-known 
models with magnetic field [H] this requirement leads automatically to the absence of the 
magnetic field. Is it possible to self-consistently consider a non-trivial electromagnetic field in 
the FLRW background? The answer is yes. Consider such a model. 

Given the structure of Einstein's field equations, in the model of parallel electric and 
magnetic fields the following electromagnetic source governs the evolution of the gravitational 
field: 



X = \(H 12 F 12 - tf 30 P 30 ) = \(M^F l2 + P 12 P 12 - M 30 P 30 - P 30 P 30 ) 

= _I (m 3 P 3 + P 3 P 3 + P 3 E 3 + E 3 E 3 ) = 1 ^M^B^ + (B^) 2 +P (3 ^E^ + (E^ 2 ) . (114) 

When both the polarization vector P^ and magnetization vector vanish, then the quan- 
tity X(t) is non- negative. Nevertheless, when P^ and M^ 3 \ are non- vanishing, X(t) can be 
negative for some time interval, or even identically vanish for the special initial conditions. 
Note that even if X is negative, the total energy density W + X is assumed to be non-negative, 
since the energy density of the material medium, W, is positive. We consider below just the 
case with X = 0. 
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5.2 Maxwell's equations 



The quantities and M ik are considered to be the function of cosmological time only. Thus, 
it follows from the second subsystem of Maxwell equations that the tensor of electromagnetic 
field has a few constant components, and using the symmetries of the model, we can choose 
only F12 = const to be non-vanishing, since the magnetic field points along the x 3 axis. Note 
that the equations ([2]) do not fix the structure of the electric field, so far an arbitrary function 
of time. We assume for the simplicity that the current P vanishes. The first subsystem of the 
Maxwell equations 

V k H ik = -^d k {^~gW k ) = -Ma 3 H i0 ) = (115) 
yj— g a 6 at 

reduces to an identity for i = and three differential equations in ordinary derivatives for the 
three components H a0 , whose solution is 

H a0 (t) = M a0 (t) + F a0 (t) = P a (t) + E a (t) = = D«(t )^M . (ii 6 ) 

cr(£) a \t) 

Note that the last relation includes the polarization vector P a (t) and electric field E a (t) only. 
Thus, the evolution of the magnetization vector M a is just governed by the constitutive equa- 
tions. 



5.3 Constitutive equations 

We use the third version of the generalized constitutive equations (11031) and (11041) . For this 
case, we specify the tetrad vectors. For the metric (11071) the normalization - orthogonality 
conditions for the tetrad vectors yield 

= 4,, = (117) 

For such a four-velocity vector the acceleration vector DU l vanishes. Consider the simplest 
relaxation model without electric conductivity, in which the polarization is coupled to magnetic 
field by the non-vanishing magnetoelectric coefficients. The master equations for such a model 
read 

r {p) V + V = (e ll -l)E- 1 B, (118) 
r (m) M + M = 7 £+ ( — - 1] B, (119) 



c 

~a{t) 



E(t)+V(t) = D(t )(^] , (120) 



'a (t 



2 



B ^- B ^{^f) ■ BM = aW)- < 121 > 

Here we use the quantities V = P^ 3 \ E = E^\ B = B^ 3 \ D = D^ 3 \ ... etc., in which, for 
simplicity, we omit the tetrad indices. Assuming that the relaxation parameters depend on 
time according to 

T( P )=6(£) , r (m) =£ a (£) , (122) 
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P( x ) = p(t )x" £ n /fi + ri (x~ 2 - x~ £ ^ 



M{x) = M(t )x «2 + T 2 [x 



x & ) + r 3 ( x 



(Ci) £2 = const) and introducing the variable x = f^jy, we obtain that the solutions to fl 1 1 8 [) 
(1T2T1) are 

(123) 
(124) 

(125) 
(126) 
(127) 



where 



D(to)(q|-l)-7^(^ 



1 



(l-26)(£||-26 



S(to) 

r 3 = 



7 2 + 



/ 1 



K^|| 



-1 (£11-26) 



-7D(t )(l-2e 



1; f > 



76 



128) 



Thus, the function X(t) reads 

x(t)=^ ( - 2 -i ) 5(t )[M(t )-r 2 -r 3 ]+^- 4 [fi 2 (to)+^ 2 (to)+5(to)r 2 - J D(to)ri] 
+^(- 2 - £ N^)[ J B(t )r 3 + J D(to)r 1 - J D(to)^(to)]. 

Since X(t) contains nine free parameters, we can choose three of them, for instance, the initial 
data, V(to), M(to) and D(to), so that X(t) = 0. In other words, we have a model with "hidden" 
electric and magnetic fields. For instance, when 

v(t ) = r u M(t ) = r 2 , (129) 

it follows that 

V{x)=Y l x- 2 , E{x)=[D{t Q )-T 1 \ x- 2 , M(x)=T 2 x- 2 , H(x) = [B(t )+T 2 } x~ 2 . (130) 
Thus, X = when 

B 2 (t ) + D 2 (t ) + B(t )T 2 - D(t )r! = . (131) 
This condition reduces to the quadratic equation for the ratio D(t )/B(t ). 



m) (i-^) (1 -2 & ) + 2 7 §m (1 _ 6 _ &)+ 



72+ GSr 26 ) (£|r26 



0. (132) 



This equation has real roots for a wide choice of the parameters £i> £2, 7, fJ>\\ and q|. For 
instance, when 6 = £ 2 = £ and 7 = two real solutions exist when — < 2£ < e\\, i.e., when 

the relaxation parameters T( p ) and T( m ) are of the order of Hubble parameter H(t)=a/a. 

Thus, we obtained an exact solution of the Einstein-Maxwell equations, describing the 
FLRW-type model, in which there is a non-vanishing magnetic field, the magnetization, the 
electric field and the polarization of the matter, however, they are hidden, i.e., their total con- 
tribution to the stress-energy tensor of the whole system vanishes. This type of behaviour was 
discussed in [10] . There the stationary magnetic field in vacuum is non- vanishing, nevertheless, 
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the exact solution to the non-minimal Einstein-Maxwell equations demonstrates the possibility 
of isotropic FLRW-type expansion. In that case the non-minimal interaction between gravita- 
tional and electromagnetic fields inspires some kind of "non-minimal screening" and hiddens 
magnetic field from the point of view of gravitational dynamics. Here we presented the example 
of "dynamic screening" , when the polarization and magnetization of the medium compensate 
the contribution of the electromagnetic field to the total stress-energy tensor. Note that in 
the proposed model the magnetic induction, = + and the electric induction, 

£)( 3 ) = + are considered to be non- vanishing. This means, that the electric polariza- 
tion compensates (partially) the contribution of the magnetic field to the expression for X, and 
the magnetization compensates (partially) the contribution of the electric field to X because of 
the special choice of the initial data. When the electric field and electric polarization are absent, 
there exists, nevertheless, the possibility that the magnetization compensates the contribution 
of the magnetic field. It is possible due to the relationship 2yU||£ 2 = 1. The magnetic induction 
is vanishing in this case. 

When X=0 the Universe is expanding isotropically and does not feel the presence of the 
electric and magnetic fields. The Einstein equations for this case are the standard (see, (11 101) 
and (lllip with X=0 and P^=P^)=P(3)=P) ■ We will not specify the equation of state and 
discuss the solutions. Nevertheless, let us note, that the application of the extended constitu- 
tive equations to the cosmic electrodynamics have something in common with inhomogeneous 
(depending on time) equations of state, introduced in [75] - [81 J . 

6 Conclusions 

The fundamentals of covariant phenomenological electrodynamics of relativistic continuous me- 
dia were elaborated three decades ago, however, the self-consistent description of a gravitating 
polarizable - magnetizable non - stationary medium is still an open question. In this paper 
we formulated extended Einstein - Maxwell model appropriated for the non-stationary electro- 
magnetically active relativistic material medium, which has three main ingredients: 

(i) . The standard Maxwell equations, describing the electrodynamic phenomena in con- 
tinuous media in terms of induction tensor and Maxwell tensor (see (pQ) and (j2J)). 

(ii) . A new set of covariant constitutive equations, containing the polarization - mag- 
netization tensor and its first, second, etc., covariant derivatives (see (|96|) . (1971) . fl 103[) and 
(11041) ). In this sense the model can be indicated as an extended one in analogy with extended 
thermodynamics. 

(in). The Einstein equations for the gravity field, in which we introduced a new effective 
stress-energy tensor (1551) . describing the contribution of the electromagnetic field, of the elec- 
tric polarization and of the magnetization of the medium. This Einstein-Maxwell model can 
be indicated as self- consistent by two reasons. First, the corresponding extended constitutive 
equations describe the interaction of matter with electromagnetic field, resulting in the dynam- 
ics of polarization - magnetization. Secondly, the polarization and magnetization contribute to 
the total stress - energy tensor, the source for the gravitational field, via the proposed effective 
stress - energy tensor of the electromagnetic field. 

Concerning our approach, we would like to emphasize three points. 

a) This extended model needs verification. In this paper we discuss only one example of the 
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exact solution to the extended Einstein-Maxwell model, describing the FLRW-type cosmological 
dynamics with hidden electromagnetic field. We prepared also a paper related to Bianchi-I 
anisotropic cosmological model with exact solutions of the new type. We hope the extended 
Einstein-Maxwell model will be also useful in application to the theory of interaction of the 
gravitational waves with electromagnetically active material media. 

b) This model admits also a generalization to gravitating static anisotropic non-homogene- 
ous media. To develop such a model one can introduce the first, second, etc. spatial derivatives 
of the polarization-magnetization tensor into the constitutive equations in analogy to the con- 
vective derivative. 

c) We introduced the extended constitutive equations by the phenomenological way. The 
next step is to confirm this approach by the consideration of the corresponding entropy pro- 
duction scalar and by analyzing the first and the second laws of thermodynamics. We shall 
present such analyses in future papers. 
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